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Wo report a measurement of the branching fraction 'Br{B'^ Di*^ oi*^) using a data sample 
corresponding to 1.3 fb"'^ of integrated luminosity collected by the DO experiment in 2002-2006 
during Run II of the Fermilab Tevatron Collider. One oi*^ meson was partially reconstructed 
in the decay Ds a-nd the other Z)i*' meson was identified using the decay Ds ^tt 

where no attempt was made to distinguish Ds and D* states. The resulting measurement is 
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Br(_Bs ^ Dg'^Dj*') = 0.039to'.oi7(stat)^p[j|Jg(syst). This was subsequently used to estimate the 
width difference APf^ in the B°-B° system: AF^^/T, = 0.079t^;E;|^ (stat)+°;°3j(syst). 



PACS numbers: 12.15.Ff, 13.20.He, 14.40.Nd 

In the standard model (SM), mixing in the sys- 
tem is expected to produce a large decay width dif- 
ference ATg = Fl — Th between the light and heavy 
mass eigenstates with a small CP-violating phase (jjs [ll- 
New phenomena could produce a significant CP-violating 
phase leading to a reduction in the observed value of 
APs compared with the SM prediction of APj/Ps = 
0.127 ± 0.024 [2]. APf^ = ^pCP-cvcn _ ^pCP-odd 

(APs — AP^"^ I cos (f>s I) can be estimated from the 
branching fraction Br(BO ^ Di*^Di*^) This decw 

is predominantly CP-even and is related to AP^^ P, [^: 
2Br(B0 ^ Di*^Di*^) « (AP?^/P,) [1 + 0(AP,/P,)]. 

Only one measurement of Br(i?^ Ds*'!?^*-*) has pre- 
viously been published, by the ALEPH [3] experiment 
at the CERN LEP collider from the study of correlated 
production of (j)(f> in decays. 

In this Letter we present a measurement of Br(i3° 
oi*'' oi*^) using a sample of semileptonic decays col- 
lected by the DO experiment at Fermilab in pp collisions 
at = 1.96 TeV. The data correspond to an integrated 
luminosity of approximately 1.3 fb^^. We present an 
analysis of the decay chain B'^ £'i*''Z3i*-* where one 
decays to (piir'^, the other D~ decays to Dj 02M~^: 
and where each (p meson decays to <^ ^ K~^K~. Charge 
conjugate states are implied throughout. No attempt 
was made to reconstruct the photon or tt'^ from the de- 
cay D* Dsj/iT^ and thus the state oi*^ oi*^ contains 
contributions from DsDs, D*Ds and D*D*. To reduce 
systematic effects, Br(S° Di*^Di*^) was normalized 
to the decay ui*^ fiv. 

The DO detector is described in detail elsewhere 
The detector components relevant to this analysis are 
the central tracking and muon systems. The DO central- 
tracking system consists of a silicon microstrip tracker 
(SMT) closest to the beampipe surrounded by a cen- 
tral scintillating-fiber tracker (CFT) with an outer ra- 
dius of 52 cm. Both tracking systems are located within 
a 2 T superconducting solenoidal magnet and are op- 
timized for tracking and vertexing for pseudorapidities 
\r]\ < 3 (SMT) and |?7| < 2.5 (CFT), where r] = 
— ln[tan(0/2)], and 9 is the polar angle with respect to 
the beam axis. The muon system is located outside of the 
liquid- argon/uranium calorimetry system and has pseu- 
dorapidity coverage \t]\ < 2. It consists of a layer of track- 
ing detectors and trigger scintillation counters in front of 
a 1.8 T iron toroid, followed by two similar layers out- 
side of the toroid. The trigger system identifies events 
of interest in a high-luminosity environment based on 
muon identification, charged tracking, and vertexing. No 



explicit trigger requirement was applied, however most 
events satisfied inclusive single-muon triggers. 

The measurement began with reconstruction of the de- 
cay chain Ds — > ^itt, 0i — s- K~^K~ ^ with tracks originat- 
ing from the same pp collision point (primary vertex) as 
a muon. All charged tracks used in the analysis were 
required to have at least two hits in both the SMT and 
CFT. Muons were required to have transverse momen- 
tum pt > 2 GeV/c, total momentum p > 3 GeV/c, and 
to have measurements in at least two layers of the muon 
system. Two oppositely charged particles with pT > 0.8 
GeV/c were selected from the remaining particles in the 
event and were assigned the mass of a kaon. An invariant 
mass of 1.01 < M{K+K-) < 1.03 GeV/c^ was required, 
to be consistent with the mass of a meson. Each pair of 
kaons satisfying these criteria was combined with a third 
particle with pT > 1.0 GeV/c, which was assigned the 
mass of a pion. The three tracks were required to form a 
Ds vertex using the algorithm described in Ref. Q. The 
cosine of the angle between the Dg momentum and the 
direction from the primary vertex to the vertex was 
required to be greater than 0.9. The Ds vertex was re- 
quired to have a displacement from the primary vertex 
in the plane perpendicular to the beam with at least 4(t 
significance. The helicity angle x is defined as the angle 
between the momenta of the Ds and a K meson in the 
(K'^K^) center of mass system. The decay of Ds <j)Tr 
follows a cos^ X distribution, while for background cosx 
is expected to be flat. Therefore, to enhance the signal, 
the criterion | cos^j > 0.35 was applied. The muon and 
pion were required to have opposite charge. The events 
passing these selections, referred to as the preselection 
sample, were used to produce the samples of (/i(/)2-Ds) 
and the normalizing sample {fJ.Ds) defined below. 

To construct a {fiDs) candidate from the preselec- 
tion sample, the Ds candidate and the muon were re- 
quired to originate from a common Bs vertex. The 
mass of the {fiDg) system was required to be less than 
5.2 GeV/c^. The number of tracks near the Bg meson 
tends to be small, thus to reduce the background from 
combinatorics, an isolation criterion was applied. The 
isolation is defined as the sum of the momenta of the 
tracks used to reconstruct the signal divided by the total 
momentum of tracks contained within a cone of radius 
ATZ = a/At/^ + Ac^ = 0.5 centered on the direction of 
the B° candidate. We required the isolation to exceed 
0.6. To suppress background, the visible proper decay 
length, defined as M{B^) ■ [Lt ■Pt)/'p't^ required to 
exceed 150 /im. Here Lt is the displacement from the 
primary vertex to the B^ decay vertex in the transverse 
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FIG. 1; (a) The {K^ K tt) invariant mass spectrum in the 
mass window 1.01 < M{K+K-) < 1.03 GeV/c^ (b) Mass 
spectrum of the {K^K~) system in the mass window 1.92 < 
M{K+K-n) < 2.00 GeV/c^ 



plane, and M{B'^) is the mass of the meson :7]. These 
data are referred to as the (fJ-Dg) sample; the result- 
ing mass spectrum of the {K^ K~ tt) system is shown in 
Fig. mi a), where the Dg and mass peaks are described 
by single Gaussians with a second-order polynomial to 
parameterize the background. The signals of Dg </>i7r 
and (/'itt"'" are clearly seen. Figure [TJb) shows the 

mass spectrum of the (K^K^) system, where a double 
Gaussian describes the (p mass peak, and a second-order 
polynomial is used to parameterize the background. 

To construct a (/i(/)2-Ds) candidate from the prese- 
lection sample, a second 02 meson from Dg (t>2l^i^ 
was required. The selection criteria to reconstruct the 
second 02 meson were identical to those of the first 

01 meson, with the exception that a wider mass range 
0.99 < M{K+K-) < 1.07 GeV/c^ was used to estimate 
the background distribution under the 02 meson. This 

02 meson and muon were required to form a Dg vertex. 
To suppress background, the mass of the (/x02) system 
was required to be 1.2 < M{^i(j)2) < 1.85 GeV/c^. The 
Ds(0i7r) and Z?s(02/^) mesons were required to form a 

vertex. The mass of the (/i02-Ds) system, i.e., the 
combined Dg — + 02 /iJ^ and Dg — > 0i7r candidates, was 
required to be 4.3 < M{p(j)2Dg) < 5.2 GeV/c^. An isola- 
tion value exceeding 0.6 and visible proper decay length 
greater than 150 /im were required for the Bg meson. 

To reduce the effect of systematic uncertain- 
ties, we calculated the ratio R — Br(_B" 



0/ii/)/Br(B^ ^ fii^Dl*'). We 



Di*^Di*^) ■ Br{Dg 
extracted Bt{B° Di*^D'g*^) from R using the known 
values J] for Bi{Dg 0/ii/), Br(B^ -> d\*'^ ^jlu), and 
Bi{Dg — > 07r). R can be expressed in terms of experi- 
mental observables: 



R 



N, 



bkg 



N, 



Di*^fii.) 



2Br 



K+K- 



eiBO 

where -/V^^Ds is the number of (fJ-Dg) events, N^^^j^^ 
is the number of (/i02-Ds) events, iVbkg is the number 
of background events in the (/i02-Ds) sample that are 
not produced by B^g decays, and /(-B^ 

D^*^ ^,v) is the fraction of events in {^iDg) coming from 
B^g D'f^fiiyX. The ratio of efficiencies £(B° 
Di*^Di*^)/e{B^ d'i ^w) to reconstruct the two pro- 
cesses was determined from simulation. All processes 
involving 6 hadrons were simulated with EvtGen Q 
interfaced to pythia followed by full modeling of 
the detector response with GEANT [10| and event re- 
construction as in data. The number of (/i-Ds) events 
was estimated from a binned fit to the {K^K^ir) mass 
distribution shown in Fig. [Ija) from the 145,000 can- 
didates passing the selection criteria. The resulting fit 
is superimposed in Fig. [TJa) as a solid line and gives 
N^,D^ = 17670 ± 230 (stat) events. 

The number of {fi(j)2Ds) events was extracted using a 
unbinned log-likelihood fit to the two-dimensional dis- 
tribution of the invariant masses Md of the (0i7r) sys- 
tem and Mff,^ of the two additional kaons from the (02m) 
system. All candidates from the (/i02£'s) sample with 
1.7 < Md < 2.3 GeV/c2 and 0.99 < M^^ < 1.07 GeV/c^ 
were included in the fit. In the fit, the masses and widths 
for both Dg and signals were fixed to the values ex- 
tracted from a fit to the (fiDg) data sample. Extracted 
from the fit were the numbers of: N^^^d^ events from cor- 
related (joint) signal production of (0i7r) and 02, events 
with a reconstructed (0i7r) in the mass peak of Z3s (0i7r) 
without joint production of 02 from (02m) (i-S., uncorre- 
lated), events with a reconstructed 02 from (02m) without 
joint production of (0i7r) in the mass peak of the Z?s(0i7r) 
(i.e., also uncorrelated) , and combinatorial background. 

The results of the fit are displayed in Fig. [2l Fig- 
ure[2Ja) displays the invariant mass distribution of (0i7r) 
candidates from the invariant mass signal window of 
£'s(02/i), and Fig. [^b) displays the 02 meson from 
Dg{(j)2^J') in the invariant mass signal decay window of 



Di*'^Di*^ 



(1) 



£'s(0i7r) candidates. The fit gives iV, 



= 13.4 



+6.6 
-6.0 



events from the 340 candidates included in the fit. 

As a consistency check, a similar sample was produced, 
but requiring the same charge for the muon and pion. 
From the fit, the number of (/i+02£'^) signal events is 
found to be zero with an upper limit of 2.6 events (68% 
CL). 

To extract the number of Bg Di*^ ^i/ and Bg 
Dg Dg events, the composition of the selected sam- 
ples must be determined. The decays Bg D^*^ ^vX 
and B°g Di*\{-^ ^jlv)vX were considered as sig- 
nal. The branching fractions for B DgD'^*^X and 
Di*^Di*'> are taken from the Ref. Q. There is 
no experimental information for the Br(B^ DgDX), 
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FIG. 2: Invariant mass distributions of (a) Ds{<j}in) events 
in the signal window of the invariant mass {K^K~) from 
Ds{4>2fj), and (b) {K^K~) events from Ds{(f2f^) in the in- 
variant mass window of the Ds{(j)in) signal region. The solid 
curve is the projected result of the unbinned log-likelihood fit, 
the dotted curve shows the polynomial background contribu- 
tion, and the dashed line shows the uncorrelated production 
of (a) Ds{4>iit) and (b) 02 mesons. 



therefore we used the value 15.4% provided by Ref. [s*] 
with an assigned uncertainty of 100%. 

In addition, the (/xDs) sample includes the processes 
cc Ds^ivX, bb — > DsfiiyX, and events with a misiden- 
tified muon, etc. which we refer to as the "peaking back- 
ground." The estimated contribution of these processes 
in the {^Dg) signal is (2 ± 1)%. In total, we estimate 
that the fraction of events in the {fJ,Ds) signal coming 



from B° DV'fii^X is /(S^ DT'fxiy) = 0.82 ± 0.05. 

We considered the number of events N^j^^^d^ from the 
(/i02^s) sample to contain contributions from 1) the 
main signal 



oi*'' Di*\ and the following back- 



ground processes 2) B -> D^*^D^*^KX, 3) 5° 



D^*'>D^*'^X, 4) -> Df'(t)fny, 5) peaking background, 
and 6) — > D^*^ ^ly combined with a <f> meson from 
fragmentation. There is no experimental information for 
most of the processes, therefore their contributions were 
estimated by counting events in different regions of the 
(/i02^s) phase space and comparing the obtained num- 
bers with the expected mass distribution for each back- 
ground process. 

The mass of the {^(t)2Ds) system for the second and 
third processes is much less than that for the main 
-)(*) n(*) 



i(*) 



least two pious due to isospin considerations. At least 
two gluons are required to produce this state (similar 
to ip{2S) J/V'tttt); it is therefore additionally sup- 
pressed and its contribution was neglected. Simulation 
shows that for the B ^ Dl Ds KX decay, the fraction 
of events with M{n(j)2D,) > 4.3 GeV/c^ is 0.05. Re- 
quiring M(/x(^2-Ds) < 4.3 GeV/c^ and keeping all other 
selections, we observe 2.8^2^g^ events in data. Assum- 
ing that all these events are due to B Dg Dl ' KX, 
we estimate their contribution to the signal {fi(f>2Ds) as 
0.14t°;^^ events. 

The fourth process produces a high mass for both the 
(^02) and (/x02-Ds) systems and requiring M(/x02) < 1.85 
GeV/c^ strongly suppresses it. Simulation shows that 
for this process, the fraction of events with M{^(f)2) < 
1.85 GeV/c^ is 0.14. Requiring M{^l(f>2) > 1.85 GeV/c^ 
and keeping all other selections, we observe 13±11 events. 
Assuming that all these events are due to the fourth 
background process, we estimate its contribution to the 
(/i02£'s) signal as 1.88 ± 1.51 events. 

The contribution of the peaking background is strongly 
suppressed by the event selection with an upper limit of 
0.4 events. We therefore included it as an additional 
uncertainty in the number of background events. 

The fitting procedure accounts for the possible back- 
ground contribution of the decay i?° ~* d\*^ to- 
gether with the uncorrelated production of a </> meson 
from fragmentation. In addition, an attempt was made 
to reconstruct {ii(f)2Ds) events in the B^ — > oi*'' fiv 
simulation containing approximately 9200 reconstructed 
(pDs) events, and no such events were found. There- 
fore the contribution of this process was neglected. In 
total, we estimate the number of background events as 
iVbkg = 2.0 ± 1.6. 

In determination of efficiencies, the final states in 
the (fJ-Dg) and (/i02^s) samples differ only by the two 
kaons from the additional 02 meson. With the excep- 
tion of the isolation criterion, all other applied selec- 
tions are the same, so many detector-related system- 
atic uncertainties cancel. The muon pT spectrum in 
dI*^ fj,i> decay differs between data and simula- 
tion due to trigger effects and the uncertainties in B 
meson production in simulation. To correct for this dif- 
ference, weighting functions were applied to all Monte 
Carlo events. They were obtained from the ratio of 
simulated and data events for pT distributions of the 
B2 meson and muon. With this correction, the ratio 



decay B^ — + Dg'Ds' because of the additional par- of efficiencies is e{B'^ 



Di*^Di*^)/e{BO, 



tides, and the requirement M(^02-Ds) > 4.3 GeV/c^ 
strongly suppresses them. The contribution of B^ —> 
Di*^Di*^X is much less than B Di*^Di*^KX because 
of higher production rates of _B+ and i?° compared to 
Compared to the B Di*'^Di*^KX process, the 
final state in the decay i?° 



Di*''D^*^X includes at 



0.055±0.001 (stat). The systematic uncertainty in this 
value is discussed below. The difference in efficiency 
is mainly due to the softer momentum spectrum of the 
muon from the Dg (p^iv decay in the (/^02^s*^) sam- 
ple, compared to the muon from the B^ — s- D^s*^ decay 
in the (/iZ?s*^) sample. 
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Using all these inputs and taking the value Br(0 
K+K-) = 0.492 ± 0.006 [3, we obtain R = 0.015 ± 
0.007 (stat). The statistical uncertainty shown includes 
only the uncertainty in N^^^d^. All other uncertainties 
are included in the systematics. 

The experimental extraction of both Br(i?g oi*^ fiiy) 
and Br(_D<, — > (j^fj-v) depend on Br(_D<, — > 4>n). Factoriz- 
ing the dependence on Br(Z?s — > (pn), we obtain Br(_B" 
^i:yDi*^)BT{D, 07r) = (2.84 ± 0.49) x lO^^^BviD, 

= (0.55 ± 0.04) • Br(i:>, -> </>7r). Using these 
numbers, we finally obtain Br(B^ D''*^Di*^) = 
0.039 j:|^:J5i^(stat). 

The systematic uncertainties in the measured value of 
Br(BO Di*^Di*^) were estimated as follows. All exter- 
nal branching fractions Q were varied within one stan- 
dard deviation. A 100% uncertainty in the number of 
background events iVbkg in the (/J.02^s) sample was as- 
sumed. The ratio of efficiencies can be affected by the 
uncertainties of reconstruction of two additional charged 
particles from the (f> meson decay. The efficiency to re- 
construct a charged pion from the decay D*^ Z?°7r+ 
was measured in Ref. ll| , and the obtained value was in 
a good agreement with the MC estimate. This compari- 
son is valid within the uncertainty of branching fractions 
of different B meson semileptonic decays, which is about 
7%. Therefore we conservatively assigned a 14% system- 
atic uncertainty (7% for each charged particle, 100% cor- 
related) to the ratio of efficiencies and propagated it to 
the final result. For the ratio of efficiencies, a 15% uncer- 
tainty was assigned for the reweighting procedure, which 
reflects the difference in efficiency between weighted and 
unweighted estimates. The dependence of the number of 
(^4>2Ds) events on the fitting procedure was estimated 
by adding a possible signal contribution from D+ events 
which decreased the correlated signal by 3%, which we 
assigned as a systematic uncertainty. 

Using these numbers, we obtain Br(i?° 



[0.039l:ogi?(stat) ± 0.014(syst)] • [0.044/Br(i:»s - . 
Using Brps •Z'tt) = 0.044 ± 0.006 [tI, we find 

Br(i?," ^ i?(*)l?(*)) ^ 0.039t;;:gi?(stat)+°:;:;^(syst). (2) 



The result is consistent with, and more precise than the 
ALEPH measurement Br(BO o'i^ D^*^) = 0.077 ± 
[4] , where the value has been recalculated us- 



0.034 



+0.038 
-0.026 



ing the current value of Br(_Ds (jiir) We calculate 
AFf^ [l| assuming that the decay B° D^*^ o'f^ is 
mainly CP-even and gives the primary contribution to 
the width difference between the CP-even and CP-odd 
_B" states 



Assuming CP-violation in 5° mixing is small this 
estimate is in good agreement with the SM prediction 
AF^/Fs = 0.127 ± 0.024 [2] and with the direct mea- 
surement of this parameter by the DO experiment in 
B° J/ipij) decays [12]. The agreement with the CDF 
measurement of AFg/Fj, which was also performed in 
J/i^P [i3|> is not as good, although still within 
two standard deviations. 
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= 0.079 



+0.038 
-0.035 



(stat)+o:o3o(syst). 
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